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ABSTRACT 


NAt  the  requaat  of  tha  Height  of  Burat  Panel,  Project  ».l  J  waa  organlaed  to  make  meaaure- 
manta  on  King  Shot,  Oparatlon  In.  that  would  aatabllah  tha  peak  ahoch  ovarpreaaure  in  tha 
blaat  vara  aa  a  function  of  dlatanca  from  tha  burat  la  tha  fraa-alr  ragion.  This  Information 
waa  required  In  particular  to  dattrmlaa  whathar  acailng  lawa  could  ba  uaad  with  axlitlng  data 
obtained  on  Oparatlon  Tumblar  to  pradlct  fraa-alr  praasuraa  from  much  tartar  waapona.  Sac- 
ondan  objective*  wara  to  racord  and  datarmlna  tha  magnitude  of  a  pracuraor  wav*  or  othar 
vlalbly  obe*r*tbU  tharmal  affacta  that  might  occur  and  to  cotlact  any  additional  Information 
that  might  explain  tha  dapartura  of  tha  fraa-alr  blaat  maaauramanta  obtalnad  on  Oparatlon 
Oraanhouaa  from  tha  Oparatlon  Tumblar  corepoclta  fraa-alr  pm* aura  raaulta.  (Tha  four  traa- 
alr  praaaura -dlatanca  curves  obtalnad  on  Oparatlon  Tumblar  acalad  vary  wall  over  tha  antlra 
pmaaum  range  measured.  Tha  compoalta  mault  la  conaldarad  to  ba  highly  rallabla.)/ 

1  Data  warn  coif  acted  uatng  the  photo -optical  technique  more  frequently  mfarrad  td 
■‘rocket  amoka-trail  photography.”  ThU  technique  haa  bean  built  around  the  ehock-ym 
method  of  paak-blaat-pr»««ura  determination.  \ 

-  The  maulu  obtained  can  taw  aummarteed  aa  followa: 

For  the  ftrabell  region  the  equation  which  waa  fitted  to  tha  radlua-tlme  data  la 

R  (ft)  -  SSOl.St*"*  (ate)  R  s  900  ft 
Tha  a  hock  overpressure  la  thla  region  la  mlatad  to  dlatanca  by 
P  (pal)  *  R'4-**  (ft)  R  a  #00  ft 


For  the  fraa-alr  region  tha  equation  that  waa  fitted  to  tha  arrival -time  data  by  tha  method 

of  leant  equarea  la 


t  (aae)  •• 


for  value*  of  tha  radial  dlatanca  R  betw 
U  In  feat  par  eecond  la 


V  •  vn.% 


i'-m‘ 


k  velocity 


#00  ft  %  R  *  3330  ft 


Tha  radiochemical  kllotoneaga  equivalent,  baaed  on  a  comp  art  eon  of  the  King  Shot  raaulta  with 
tha  Operation  Tumbler  eompoatla  raaulta,  U  Ml  *  30  kt  Tha  TNT  efficiency  of  King  Shot  waa 
found  to  be  II  I  •  3  0  par  seat  within  tha  preeaure  range  of  300  to  30  peu 


xThe  precursor  were  detected  orsr  the  land  im  bp  pressure -Oms  gauges  of  Project  6.1 
(WT  -601)  was  not  observed  to  the  mottos  -picture  (lias  exposed  particularly  (or  Protect  (.13. 
Palm  trees  on  an  Island  la  Um  foreground  obscured  the  Island  nearest  the  burst,  where  the  pre¬ 
cursor  wave  was  detected.  No  precursor  was  observed  over  the  water  either  bp  gauges  or 
photographically.  No  positive  evidence  was  (ouod  to  explain  the  departure  o(  Operation  Green¬ 
house  free -air  measurements  from  the  Tumbler  composite  tree-air  pressure  results.-  .- 
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CHAPTER  1 


INTRODUCTION 


1.1  OBJECTIVES 

Upon  the  recommendation  of  the  Height  o f  Buret  Panel,  the  U.  I.  Naval  Ordnance  Lahore* 
lory  (NOL/  waa  requested  by  the  Armed  Forcee  Special  Veepocs  Project  (AFSWP)  to  parttcl- 
pate  In  Eng  Shot,  Operation  Ivy.  In  particular,  the  NOL  van  asked  to  conduct  the  neeeasary 
experiments  required  to  ohtein  the  following  information: 

1.  Peak  shock  overpressure  ae  e  function  of  distance  in  the  free-air  legion.  This  Infor- 
mail  on  was  required  ia  particular  to  determine  whether  seeling  lawe  could  be  used  with  exist¬ 
ing  date  obtained  on  Operation  Tumbler  to  predict  free-slr  p raeeuree  from  much  larger 
weapons. 

S.  Information  relative  to  the  formation  sad  magnitude  of  any  precursor  wavs  or  other 
visibly  obeerreblf  thermal  '•fleet. 

I.  Additional  information  that  might  explain  the  departure  of  the  free-air  measurement* 
obtained  on  Operation  Qrevnhoue*  from  the  Operation  Tumbler  compoeite  free-air  pressure 
resuUa.  (The  four  free-air  preaaure-dieUnce  rurvre  obtained  on  Operation  Tumbler  sceleo 
very  well  over  the  entire  preeeure  range  measured.  The  compoeite  result  it  considered  to  be 
highly  rellabie.) 

1.1  BISTORT 

b  order  to  cbtatn  the  detlred  Information  required  by  objective  I,  the  photo-op<iest  tech- 
tugwe,  which  ta  tow  Wwa  ae  the  "rocket  emoke-U-ait  photography  method,"  seemed  to  be  ade¬ 
quate  with  only  e  alight  modification.  This  technique  wie  developed  at  the  NOL  for  ua*  la 
Operation  Oreenhouee  where  It  proved  to  be  quit*  cucreie'vl.1  Highly  eetin/ectory  result* 
were  also  obtained  later  os  Operations  Jangle*  and  Tumbler  1 

On  the  fourth  shot  of  the  Tumbler  oenee  of  leetr.  a  ?r .cursor  wavs  appeared  is  ihe  roeset 
smoke-trail  films  sad  was  naatyeed  futAtUellvely  for  the  first  time  (Suck  waves  si»o  eg 
peered  la  certata  films  of  Operation  Butter,  ta  which  '.be  NOL  did  not  participate.  but  were  sot 
analysed  at  the  time.)  Thus  tt  wts  believed  the',  l«  the  tarn*  aig.trlmeni  designed  principally 
to  aceompliah  otjeciiv*  I,  the  information  sough*  objective  I  would  be  obtained 

Wiih  regard  to  the  test  objective  the  rock*',  moke -trait  method  -eed  in  the  four  eh  a*  «f 
Operat.ua  Tumbler  yielded  euck  conels’eit  measurement*  after  the  eppllraUon  o«*  ere  ic*  lee* 
that  a  highly  reliable  p re emire -distance  curve  vie  rad*  eve  table  Comp*/ lent*  uf  Tumaiei 
and  Oreenhouee  frat-alr  dele*  showed  that  the  prteeuree  ebeerved  on  Greenhouse  were  in 
creasing  ?  hightr  then  thoee  obtained  on  Tumbler  *•  ike  equivalent  reduced  dleiev*  frrvw  ike 
brweh  tacreeaed  P  reeled  check  ig  of  all  Ihe  date  ne.  tied  Ike  poeeibillly  of  error  l»  ike 
Oreenhouee  celculatloee,  and  It  vte  shown  con*- lualmiy  that  tke  arrival  ilniee  observed 


\ 


Operation  Greenhouse  led  to  higher  shock  velocities,  t nd  hence  higher  pressures,  et  equivalent 
seeled  dletencee 

On  Operation  Ivy  it  wee  hoped  that  some  Information  might  be  gstned  to  solve  this  problem. 
Teste  on  Opert'lon  Tumbler  vers  conducted  et  the  AZC  Nevada  Proving  Grounds  (NPO)  at  en 
ettitude  o t  approximately  4000  ft  abere  see  level,  whereas  Operetton  Greenhouse  took  piece  el 
Enlwetok  Atoll  In  the  Marehell  lalenda,  virtually  et  aee  level.  The  etmoephere  et  the  NPG  wee 
dry,  whereas  the  lelend  site  wee  comperetlvely  motet  and  humid.  It  wee  thought  thet  possibly 
the  cause  of  the  disagreement  might  stem  from  these  atmospheric  differences  or  from  the  feet 
thet  the  Greenhouse  teste  were  (ov  er  shots  wherees  thoae  of  Tumuler  were  elrdrope 

During  the  planning  phese  It  wee  ennounced  that  King  Shot,  Operation  Ivy,  wee  to  be  an 
airdrop  similar  to  thoae  on  Tumbler  but  et  the  Island  site.  This  presented  an  excellent  oppor¬ 
tunity  to  resolve  'he  difficulty  at  least  partially,  If  Indeed  such  environmental  conditions  were 
the  cause  of  the  disagreement. 

IS  OPERATIONS 

Project  6. IS  wee  organised  In  two  groupe,  one  which  went  Into  the  field  and  one  which 
aided  In  the  preteet  preparation#  and  poetteet  analytic  et  the  NOL.  Administrative  details 
within  the  project  were  carried  out  jointly  by  J.  t  Moulton,  it.,  and  P.  Hanlon.  E.  P.  Cox  and 
F.  B  Pores!  were  the  directors  of  Program  6  and  provided  over-ell  supervision  and  technical 
guidance. 

Project  8. IS  pereonncl  Included  i.  T.  Moulton,  Jr.,  Project  Officer,  Aneiyrle;  the  field 
party  composed  nf  P.  Hanlon,  field  Project  Officer,  Analysis;  B.  M.  Lorlng,  Supply  and  Anel- 
yele;  end  C.  L.  Karmel,  Inetrumentatlon  and  Analysis;  end  J.  R.  Mitchell,  Supply  Officer.  In 
the  field,  Karmel,  as  the  assistant  to  the  field  Project  Officer,  was  largely  responsible  for 
Instrumentation .  Lorlng  served  as  Sipply  Oiflcer  In  lhe  field  and  esaleted  lr.  the  lnatellatlon 
of  equipment.  Labor  wee  provided  by  Holmte  A  Nerver  through  Teak  Group  (TO;  133.1. 

1.8  BACKGROUND  DETERMINATION  OP  PEAK  SHOT  OVERPRESSURE  BY  THE 

VELOCmr  METHOD 

The  photo-optlcel  technique  uaed  on  King  Shot,  Operetton  Ivy,  to  obtain  arrival-time  date 
for  the  determination  of  peak  shock  overpressures  le  esqilaJned  In  detail  In  reporta  on  Opera¬ 
tion#  Greenhouse1  and  fugle. 1  Briefly  the  technique  eoneleted  In  eetebliehing  e  rocket  emoke- 
trell  grid  behind  the  buret  and  recording  'he  ahock-wevt  growth  ee  e  function  of  time  with 
high-speed  motion-picture  cameras  The  position  of  the  shock  front  te  determined  by  record¬ 
ing  photographically  the  light  reye  reflected  from  Lhe  grid;  thoae  rays  which  peas  tengentlelly 
to  the  shock  front  are  refracted,  causing  bretha  tc  appear  to  the  otherwlee  continuous  grid 
Unee.  With  Lhe  explosion  certs;  as  aero,  frerae-by -frame  measurements  of  the  distance  ere 
made  on  30  *  magnified  Images  In  e  direct  -projection  Rcccrdck.  Time  per  frame  le  also  noted. 

Prom  these  errlvel-tlme  data  the  instantaneous  shock  velocity  et  known  dletencee  can  be 
determined  This  la  beet  dors  by  fitting  the  data  with  e  smooth  curve  which  car.  be  expressed 
In  rloaed  mathematics!  form.  Differentia? 'on  of  the  equation  yields  an  expression  for  the  ve¬ 
locity  as  e  function  of  radial  distance  from  the  buret. 

In  the  early  stages,  l.e  ,  the  fireball  region,  an  exponential  function  le  uaed 

R  -  Kl*  (1  l) 

whc.  H  •  distance  from  weapon  aero 
t  ■  time 
K  *  constant 

n  ■  the  elope  of  the  log-log  plot  of  R  and  t 
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l>tfferentisli<m  erf  lhl i  equation  fives 
D  - «W*  1 


(1  >' 


where  t  11  ihe  ihetanlaneoui  shock  velocity 

Or  prior  ttsti  (h*  free-nr  arrival-time  dm  were  fitted  by  i  cubic  polynomial  by  th« 
method  of  lull  squares  On  this  t»  new  fitting  function  It  Introduced  The  new  equation  li 


b 


rR-TldRtC 


(1.3) 


where  «,  b,  end  C  ere  constants.  At  first  (line*  thli  equetlon  appeera  dumiy,  but  It  la  readily 
fitted  by  the  method  of  leeet  aqueree  to  the  errtvel-tlme  dale  on  IBM  equipment.  The  advan¬ 
tages  of  the  uae  of  the  new  equation  are  twofold:  (1)  A  more  realiatlc  approech  to  the  eolutlon 
of  the  problem  la  made  In  view  of  the  physlcel  nature  of  the  phenomenon,  end  (2)  conalderable 
time  le  saved  in  the  analyili  of  thie  type  of  data.  Upon  dlfftrenilallon,  e  comparatively  simple 
function  results: 


“■•MU"]  11,1 

Fitting  the  King  Shot  data  by  both  the  new  form  and  the  polynomial  ueed  previously  led  to 
Identical  pressure  results.  (A  comparison  of  the  Initintineoue  shock  velocities,  which  la  an 
extremely  critical  comparison,  Indicated  that  agreement  between  the  two  nets  of  results  was 
T'lthln  *0  34  per  cent  over  the  entire  range.)  Some  of  the  Tumbler  data  have  been  fitted  using 
the  new  form,  and  a  comparison  showed  the  pressure  results  of  both  methods  to  be  Identical. 

A  full  explanation  of  the  derivation  of  the  equation  and  the  method  of  fitting  date  by  the  new 
function  are  given  in  the  Appendix. 

The  peak  pressure  In  the  shock  wave  le  a  known  function  of  the  velocity  and  can  be  cal¬ 
culated  using  the  Rankins -Hugonlot  relation:* 


where  P*  ■  peak  ahock  overpraaaure  (pal) 

P,  •  ambient  pressure  ahead  of  the  ahock  (pal) 
y  •  the  ratio  of  specific  heats  for  air,  1.40 
C,  *  speed  o<  sound  ahead  of  ahock  (ft  sec) 

U  »  .natanUnvoua  ahock  velocity  (ft/eec) 

P,  and  T«  are  measured  directly,  and 


*tn  regions  of  very  high  pressure  the  Fanklne-Hujonlot  relation,  aa  written,  la  In  error 
because  the  equation  of  state  on  which  It  la  based  no  longer  applies.  Furthermore  the  ratio  of 
specific  heala  for  air  (>)  becomes  meaningless  To  overcome  the  difficulties  Introduced  by 
these  variations,  corrections  are  made  by  uae  of  the  Hlr achielder -Curttar  tables'  which  give 
P,  In  terms  of  l'.  the  changes  In  y  and  the  equation  of  aisle  being  taken  Into  account  (are  also 
Sec  3  !,• 


13 


irwmers 


I',  -  INI^I  '  ^ 


ii  •> 

•  Mra  T.  I*  the  (MtHWtl  tcmparetura  In  i*|fNl  «e«ti|re4a  end  C,  II  Ik*  *al«ily  nt  e<«und  In 
(••I  par  refund 


REFERENCE! 

I.  J  7  Moult  on,  Jr  ,  and  B  T.  RbbMi,  Peak  Praaaura  va  Diet  ant  a  la  tha  Free-air  an!  Mach 
Raglan  a  Uetng  fmofca-rocket  Photograph*.  Oreenhouae  Report,  Anne*  1  I,  pan  (I,  tar  I, 
WT  M. 

I.  J.  F.  Moulton,  Jr.,  K.  R.  Walthall,  and  P  Hanlon,  Path  Praaaura  ve  Dletance  In  Fraa  Air 
Uaing  laoka  Rocket  Photography,  Buatar>Jangla  Project  l  lb  Report,  WT  311 

I.  C.  J.  Aronaon,  J.  7.  Moulton,  Jr.,  at  at.,  Fraa  Air  and  0  route.  Laval  Praaaura  Maeeuia- 
■anta,  Tumbler -Ss  sppar  Pr<i acta  M  and  1.1  Report,  WT-Sll. 

4.  J.  O.  Hirech/aldar  attd  C.  7.  Curtlae,  Tharmojyc anile  Propartlaa  of  Air,  Vol.  □,  Cm »a rally 
of  Wiaconaln  (NHL),  Dae.  It,  IMS. 


14 


e-1* -i”.’ «'’■■■»  "»  W  ■'h!l"w,aa 

•W  J  V'^TS*  W  W  W 


i  v  ?  o . 


CHAPTER  3 


INSTRUMENTATION 


A  plan  view  of  tha  Project  6  13  Instrumentation  layout  la  shown  In  Fig  3.1. 

3.1  ROCKF  T  SMOKE-TRAIL  DETECTION  GRID 
a. 1.1  Smoke  Rockets 

To  establish  the  shock-wave  detection  grid,  18  nmoke -producing  rockets  wete  fired  7  aec 
prior  to  burst  time  from  Station  8140.  Each  round  waa  s  modified  3*0  Rocket  Head  Mark  10 
with  a  3*0  Spin  Stablllied  Motor  Mark  3.  Tan  pounds  of  FS  chemical  amoke  mix  was  re¬ 
leased  from  each  head  during  ita  upward  trajectory.  The  amoks  rocket,  as  used  on  this  and 
previous  operations,  waa  developed  for  this  purpose  at  the  NOL1  In  1830. 

a.l. a  Rocket  Launchers 

On  previous  operations  the  amoke  trails  were  established  In  the  form  of  a  vsrtkaMlns 
grid.  Thta  form  waa  not  feasible  In  ths  Ivy  teat  because  of  limited  dry  land  areas  on  which  to 
locate  rocket  isuncnsrs.  in  its  place  ■  fen-type  grid  was  used,  ss  was  first  suggested  by  C  J 
Aronaon  of  the  NOL  (see  Fig.  2.3).  The  rocke*  launchers, 1  in  two  idsntlcal  ba 1 1 cries  facing  in 
opposite  directions,  ware  aligned  ao  that  the  direction  of  fire  was  perpendicular  to  the  llns  of 
sight  of  the  cameras  (Figs,  a.l  and  3.3).  Each  bsttsry  conalsted  of  nine  launchers  with  elevation 
angles  of  10.  30,  30,  40,  50,  60,  10,  80,  and  83*.  At  one  end  of  the  launching  plot,  s  single 
launcher,  st  60*  elevation,  was  arranged  to  fire  parallel  to  the  line  of  sigh!  (Figs  3  4  end  2  3) 

3.1.3  Power  and  Timing 

Considerable  saving  In  construction  and  lsbcr  costs  was  effected  through  the  use  of  a 
novel  power-and-tlralng  station  Timing  relays,  ttms-delay  devices,  and  power  take-off 
switches  wars  mounted  Inside  a  waterprooied  wood  an  box  the  else  of  a  standard  fool  locker. 
Before  evacuation  the  necessary  switches  were  closed  In  the  ready  position,  and  Ihe  box  waa 
sealed  and  burled  beneath  sandbags  The  station  provsd  to  be  just  aa  satisfactory  as  ths  re- 
lnforced-concrete  blast  sheltsrs  coating  ssveral  thousands  of  dollars  that  were  used  In  pre¬ 
vious  operations 

The  circuitry  and  all  othar  details  concerning  ths  aatahllahment  of  (ha  rocket  smoke-trill 
grid  are  given  in  reference  3 

3  3  PHOTOGPAPHIC  INSTRUMENTATION 

The  photographic  records  (or  Project  8  13  were  obtained  by  lha  jta/f  of  Edgerlon,  Or 
meshauaen  (  Grier,  Inc  (EGtG)  Three  high  speed  Mitchell  cameras  were  Instilled  in 
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lhation  306  for  the  specific  us*  of  thi*  pro.'*ct,  b'it  on*  of  them  failed  to  t anctlon.  Comp’ete 
photog«aphlc  details  sr*  given  in  Table  3.1 

Th*  use  cf  fiducial  markers  via  unnecessary  on  this  t* at  Th*  bat*  lln*  of  the  plane  of 
meaauremen*  was  readily  established  by  utillilng  th*  horlton.  On  Operation  Tumbler*  II  was 
determine*,  that  all  that  Is  required  for  dlatarc*  calibration  la  an  accurate  calibration  "rcae” 
placed  on  tht  film  before  processing,  the  accurately  measured  focal  length  of  th*  lens  used, 
and  the  rang*  to  th*  desired  objective  plan*.  Vertical  and  horliontal  scales  can  be  established 
with  sn  sccurscy  of  better  than  0.01  per  cenl  In  this  fashion. 

Timing  msrks  were  plsced  on  th*  film  during  the  recording  period  st  th*  rale  of  1  36  cpa. 
This  very  low  rst*  is  discussed  further  In  Secs.  3.3  and  3.7.2. 


Tsbl*  2.1  -PHOTOGRAPHIC  DETAILS 


Film  No. 

Station 

Camera 

Effective  aperture 
Effective  focal 
length,  mm 
Nominal  frsm*  rate, 
frame  s/aec 
Timing  marts,  cps 
Vertical  aiming 
Horlsontal  aiming 


16283 

16261 

306 

306 

MMH-1 

MMH-7 

f/ii-is* 

f/11-15 

88.86 

35.28 

85 

80 

1  86 

1.96 

4*  OS’  4’ 08' 

Gif*  GZ» 


•Ground  aero. 
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CHAPTER  3 

RESULTS 


3.1  INSTRUMENTATION  AND  RECORDS 

The  rocket  smoke-trail  fan  grid  was  established  just  prior  to  aero  time  and  was  satis¬ 
factory  in  most  respects  As  can  be  seen  In  Etg  2  2,  the  grid  appeared  to  be  somewhat  lacking 
in  contrast  In  the  upper  right-hand  portion  of  the  object  plane,  but  this  did  not  interfere  seri¬ 
ously  with  measurements  made  on  the  film  originals  If  the  smoke  trails  had  been  approxi¬ 
mately  1000  ft  closer  to  the  burst,  the/  would  have  appeared  in  greater  contrast  but  possibly 
would  have  been  in  the  region  where  trail  evaporation  occurs  The  grid  was  placed  at  the  dis¬ 
tance  chosen  so  that.  If  the  atmosphere  were  clear  (greater  than  98  per  cent  transmission), 
good  contrast  could  be  expected. 

Of  the  three  cameras  assigned  specifically  to  Project  8  13,  only  two  functioned  properly. 
The  camera  having  a  SO-mm  focal-length  lens  failed  to  operate.  Records  were  obtained  by  the 
cameras  having  lenses  of  100-  and  35-mm  focal  length,  however,  only  the  record  from  the 
earners  equipped  with  the  100-mm  focal-length  lens  was  used  for  the  analysis  The  other  rec¬ 
ord,  which  should  have  been  of  greater  value  from  the  standpoint  of  larger  usable  field  of 
view,  proved  to  be  of  little  additional  quantitative  value  because  of  the  poor  atmospheric  con¬ 
ditions  that  existed.  The  high  humidity  led  to  such  low  contrast  that  the  rocket  smoke  trails 
were  not  visible  much  farther  away  from  the  center  of  burst  on  the  3S-mm  film  with  a  wtdi  r 
field  of  view  than  on  the  longer-focal  -  length  record  This  condition,  coupled  with  poorer  space 
resolution,  made  the  short -focal-length  record  of  qualitative  value  only 

3  2  TIME  AND  DISTANCE  SCALES 

.he  arrival-time  data  for  King  Shot  were  measured  from  film  16293  The  vertical  plane 
of  measurement  Is  shown  In  Fig  2.1  The  timing  marks  placed  on  the  film  at  the  rate  of  1  96 
cps  were  not  sufficient  by  themselves  to  ascertain  whether  the  film  speed  was  cons.ant  during 
the  period  of  Interest  However,  from  additional  Information  available  It  was  determined  that 
In  all  probability  the  film  speed  was  constant  The  time  per  frame  was  found  to  be  0  01 1 7N4  » 

0  0C0150  sec  (the  last  two  decimal  places  were  held  for  computational  reasons'  Timing  accu¬ 
racy  is  discussed  in  detail  in  Sec  3  7  2 

The  distance  scale  for  both  horizontal  and  vertical  measurements  was  found  to  be  10  686 
ft  mm  on  the  eilar,  d  image  Thus  measurements  made  to  0  1  mm  with  relative  ease  wrre 
far  better  than  ihe  maximum  sialic  resolution  uncertainly  of  >6  ft  Accuracy  considerations 
are  developed  further  In  Sec  3  7  3 

3  3  AHRIVAI  TIME  DATA 

Arrival-lime  data  were  obtained  over  the  ranges  600  to  3250  ft  In  distance  and  0  0  to  0  7 
sec  In  lime  These  data  are  guen  in  Table  3  1  and  are  plotted  in  Fig  3  1  toy  'thrr  with  the 


curve  fitted  to  th«  dett  by  the  method  a!  Itiit  tquerei.  A*  ctn  be  eeen,  the  ecetter  tn  th»  dm 
te  email 


T»bl»  SI  —  Alt  RIVAL- TIME  DATA*! 


T 


Dtatance 
from  weapon 
aaro,  ft 

Arrival 
ttma,  aac 

Dletance 
from  wetpjn 
aero,  ft 

Arrival 
time,  etc 

Dletance 
from  weapon 
aaro,  ft 

Arrival 

time,  eec 

680  SIS 

0.011320 

3013  463 

0.2464  00 

2693  618 

0  481480 

781  SIS 

0.023074 

2048.872 

0.268154 

2732  380 

0  493234 

897.800 

0.03462C 

2088.886 

0.289908 

2749.674 

0  604986 

1000.087 

0.048582 

2137.941 

0.281662 

2791.133 

0  51-743 

1094. 70S 

0.058338 

3183.091 

0.393416 

38  04  514 

0  628496 

1175  607 

0  070090 

3303.964 

0.303170 

2843  958 

0  540260 

1261.078 

0.081844 

3344.839 

0.31S924 

3866  850 

0.552004 

1320  138 

0.093  6  98 

3287.949 

0.326678 

2900  026 

0  383758 

.  •  . 

•  .  • 

3316.206 

0.340432 

3921  398 

0.5755’  3 

1469.711 

0.117106 

3145.628 

0.353183 

3943.334 

0.587266 

1538.064 

0.138860 

2385  539 

0.363940 

21*74.888 

0  599020 

1678.978 

0.140614 

2418.938 

0.373694 

3003  688 

Q.610"74 

1633.078 

0.153368 

2443. 006 

0.387448 

3037  851 

0  622528 

1880.685 

0.164113 

3486.919 

0.399203 

3058  334 

0  634282 

1744.081 

0.175878 

3513.489 

0.410958 

3093.487 

0.64  6036 

1799.049 

0.187830 

3649.803 

0.433710 

3113.091 

0  457790 

1843.263 

0.199384 

3580.963 

0.434464 

1145  188 

0  609544 

1874.989 

0.311138 

3811.163 

0.448218 

3183  084 

0  661299 

1919.371 

0  313893 

3843  579 

0.457973 

1969.378 

0.334648 

1683.018 

0.4897*4 

3313  146 
3340.398 

0  704606 

0  718560 

•All  decimal  plica*  In  dlatancen  and  tha  lut  two  decimal  plicaa  In  tlmaa  art  held  tor  com¬ 
putational  rsaeone  only. 

TTh#  taro  Urnaa  on  the  rocket  trill  film*  were  aatabllahad  by  prelim. nary  fiitball  dm 
meuured  at  the  Teel  Site.  Later  and  more  complete  meaeuremt.nte  have  re  lulled  in  radiue- 
tlzna  cum  change*  which  give  an  error  of  0  25  maec  In  tha  abeoluta  time  of  the  radlua-tlme 
Currae  ehown  bare.  i’ecaue*  of  the  ecetinj  mat  bod  a,  thle  abeoluta  error  doe*  not  affect  the 
validity  of  tha  preeewe-dletance  curve  or  tha  tealed  yield. 


13.1  Firaball  Region 

In  the  early  atagea  of  ahock  growth,  la,  the  fireball  region,  tha  radlua  increased  expo¬ 
nentially  with  time,  and  the  data  form  •  atralght  line  whan  (lotted  on  a  leg  log  graph  a*  In 
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Fig.  3  3.  The  curve  fitted  to  thtat  data  K 

R  •  3 30 3  3 1  * ,M  600  ft  s  R  s  900  ft  (3.1) 

•her*  R  l*  the  dlatanct  from  weapon  atco  In  ftet  and  t  t»  the  ttmt  In  atcond*.  Tht  flraball 
raeaaurementa  mad*  on  tht  roeket-amokt-trall  photograph*  art  In  txctllant  agreerajnt  with 
thoat  madt  by  EG40  on  tht  hlgh-apttd-camtra  rtcorda  of  conalderabty  bttttr  tlmt  rtuolu- 
tlon  Aa  dlacuaatd  In  Stc  3.7.1,  thla  fact  waa  uatd  In  tht  timing -accuracy  dttarminatlcvi 

3  3  1  Frtt-alr  Rigion 

Ovtr  tht  Iptttr  portion  of  tht  arrlval-tlmt  curve,  tht  data  wtrt  fltttd  by 

fo  oom^-j  or|  -  0.6136  (3.1) 

•htrt  R  la  tht  radial  dlatanct  from  wtapon  ttro  ovtr  tht  rangt  900  ft  s  R  s  3130  ft. 

Mtaaurtmtnta  of  tha  fir  at  ftw  framtt  of  tht  film  rtcord  showtd  that  tht  flraball  growth 
•aa  aymmttrlcal  about  tht  burat  point.  Shortly  afttr  ahock  breakaway,  a  ftw  framta  Indicated 
that  tht  ahock  wave  out  to  about  1300  ft  waa  alto  aymmttrlcal  about  the  ctnttr  of  tht  burat, 
but.  In  tht  latter  thrtt-fourtha  of  tht  rtcord,  only  tht  taattrn  aide  of  the  fan  grid  vaa  vlaiblt. 
Aa  a  rtault  It  cannot  bt  aatd  with  certainty  that  the  ahock  wave  waa  aymmttrlcal  ovtr  thla 
region. 


3  4  METEOROLOGICAL  DATA 

Tht  atmoaph  <c  prtaaurt  P,  and  ttmperaturt  T,  wtrt  mtaaurtd  prior  to  tht  ahot  Tht 
datt  from  both  au.  .act  and  upptr-alr  obatrvatlona  wtrt  uatd.  Tht  velocity  of  aound  C,  at  the 
variou*  Uvtia  waa  computtd  ualng  Eq.  1.6. 

Both  P ,  ind  C„  for  lack  of  more  complete  data,  wtrt  aaaumtd  to  vary  linearly  with  alti¬ 
tude  Thtat  datt  art  gtvan  In  Table  3.1  and  plotted  In  Fig.  3.3. 


Table  3.1 —  METEOROLOGICAL  DATA 


Altitude, 

ft 

Prtaaurt, 

pal 

Temperature, 

•c 

Velocity 
of  aound, 
ft /ate 

0 

14.649 

18.9 

1143  9 

310 

14.804 

11.8 

1143.8 

1440* 

13. 96* 

15.4* 

1138  !• 

1390 

13.833 

33.3 

1133  9 

'Interpolated  from  Radtoaonde  data. 


3  3  PEAK- SHOCK- OVERPRESSURE  -  DISTANCE  DATA 

Equ.t  lone  3  1  and  3  1  were  differentiated  to  obtain  expraeelor.a  fo;  the  Inattnlar.e out  ahock 
veloctliea,  U.  al  dealred  dlatanct*.  For  the  fireball  region 

U  -  1381  Jf*  •'*  600  ft  a  R  *  900  ft  (3  3) 
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PlrabaU  ra<flua  rt  dot,  NOt  aad  EGfcG  data. 
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Is  obtetned  This  r»»ult  leads  to  the  relation 
P  *  FT*  ,3 

which  Is  in  effective  agreement  with  theory.*  For  the  free-slr  region 
U  -  707.5  |l  ♦  800  ft  J  B  s  3250  ft 


800  ft  a  R  s  3250  ft 


These  values,  togsther  with  the  corresponding  meteorological  dits  from  Fig.  2.3,  were  used  to 
enter  the  Hirschfelder-CurtlM  tables1  from  which  the  peak  shock  overpressures  ware  obtained 
for  the  corresponding  distances.  The  values  are  given  In  Table  3.3  and  plotted  In  Fig.  3.4.  For 
pressu'es  below  approximately  100  pal  the  values  from  the  HSrschfelde'-Curtlss  tables  agree 
exactly  with  those  that  would  be  obtained  using  the  Rsnklne-tfugcnlot  pressure-velocity  rela¬ 
tion,  Eq.  1  5.  Above  the  100-psl  level,  variations  In  y  and  the  equation  of  state,  on  w.,ich  the 
Ranklne-Hugonlot  relation  depends,  begin  to  take  significant  effect.  These  variations  are  ac¬ 
counted  for  In  the  Hlrschfelder-Curtlss  tables.  For  these  reasons  values  from  these  tables 
hsve  been  used  throughout  the  calculations. 

Table  3.3  — PE  AK-OVER PRESSURE  -  DISTANCE  DATA 


A  word  about  the  use  of  the  meteorologlcsl  date  In  carrying  out  the  prssaurs  calculations  Is 
In  order.  For  the  sake  of  uniformity  with  other  air  bias'  dels1  published  previously,  the  values 
of  P,  and  C,  were  taktn  slong  a  vertical  line  from  weapon  aero  to  ground  .ero  (CZ)  at  dis¬ 
tances  corresponding  to  those  selected  for  substitution  In  the  shock-velocity  equations.  Tor  all 
distance*  greeter  then  the  burst  height  (1480  ft),  the  values  ooaerved  st  an  altitude  of  100  ft 
wc  *  used. 

3  8  PRECURSOR  WAVE  AND  THERMAL  EFFECTS 

The  excellrnt  prseeure-tlms  records  obtained  by  Project  6  1  (*ee  WT-60J)  Indicated  con¬ 
clusively  the  existence  *  precursor  wave  over  the  'vnd  eras,  whereas  (eugee  pieced  over 
the  water  detected  a  standard  pteture-book-type  shock  wave  Palm  trees  on  an  lslend  In  the 
foreground  of  the  am  xa-trall  photograph*  complately  obacurad  GZ  and  th»  Itland  whare  the 
precursor  wee  detected  No  precursor  vn  observed  over  the  water  eurfece,  which  wee  clearly 

•Theorlee  posed  by  G  1  Teylor,  J.  G  Kirkwood  »nd  S  R  Brinkley.  T  B  Porte!,  end 
others  ere  In  essential  egreetr.en'  with  regard  to  the  value  of  the  exponent  In  thte  reietion 
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visible  in  the  films  This  serves  to  confirm  ths  gsuge  results  over  wster.  No  thermal  effects 
of  any  kind  wsro  observed  In  these  photographs. 

3  7  ACCURACY  OF  RESULTS 

3.7.1  Sources  of  Error 

The  possible  sources  of  srror  and  the  procedures  for  calculating  their  magnitudes  are 
discussed  In  detail  In  reference  3  Each  of  the  seven  major  sources  of  error  listed  below  were 
given  due  consideration: 

1.  Timing  calibration. 

2  Static*  and  dynamic -resolution  uncertainties  associated  with  film  measurements  under 
Ideal  conditions. 

3.  Scaling  distance  on  film. 

4.  Forsshortsnlng  of  the  Image  In  the  plane  of  meseurement. 

5.  .Meteorological  data. 

S.  Curvs  fitting. 

7  The  variation  of  y. 

3.7.2  Timing  Accuracy 

The  small  number  of  timing  marks  which  wars  placed  on  the  film  at  a  rate  of  1  96  cps 
mads  It  difficult  to  establish  quantitatively  ths  uncertainty  In  'hi  time  resolution  ft  was  Im¬ 
possible  to  determlns  directly  whether  ths  film  speed  was  constant  over  the  entire  region  of 
tntsrsst.  The  average  film  speed  was  measured  whersvtr  possible.  In  ths  region  of  Interest 
and  the  regions  lmmdlstely  beyond,  u.'  film  speed  wss  found  to  bs  the  same.  The  first  six 
frames  of  the  record  (film  1C283)  Included  the  fireball  and  early  transition  regions.  In  which 
exceptionally  accurate  radial  distances  could  be  measured.  Vhen  the::  wars  compared  with 
ths  fireball  growth  measurements  of  EG&G,  for  which  timing  wss  highly  sccursts,  It  was  found 
that  ths  agreement  was  excellent.  Thle  la  shown  In  Fig.  3  2.  Ae  •  retuu  the  time  per  frame 
was  eatabllehed  accurataly  over  this  region  to  within  *0.000130  aec  per  frame.  It  wea  assumed 
that  ths  film  rate  wee  constant  and  that  the  time  per  frama  uaed  in  the  fireball  region  waa 
valid  throughout  ths  entire  period  of  Interest.  In  the  subsequent  enslysis  the  resulting  data 
Indicated  that  the  film  speed  wea  constant  by  the  more  or  less  random  distribution  about  the 
fitted  arrival-time  curvs  of  Fig.  3.1.  In  addition,  If  ths  momentum  of  the  film  end  the  film 
drlvs  is  considered,  It  seems  unlikely  that  there  would  be  ■  significant  vsr.atton  in  speed  with¬ 
in  0.7  sec.  The  assumption  of  constant  frame  rets  le  thus  considered  reasonable. 

3.7.3  Distance  Accuracy 

The  accuracy  In  scaling  distance  la  dependent  on  the  measurement*  of  the  focel  length 
cf  the  lens  and  the  distance  from  the  earner*  to  the  object  plan*.  The  uncertainty  In  these 
rr-aaureme  ts  1*  known  to  be  less  than  0.1  per  cent.  After  the  center  of  burst  le  located  on  the 
film  In  any  given  vsrly  frame,  its  position  Is  determined  In  the  later  frames  by  fixing  its  posi¬ 
tion  with  rasped  to  the  film  sprocket  holes.  Assuming  tl. '  earners  to  be  operating  normel’y, 
the  variation  of  the  position  of  tht  sprocket  hole*  with  respect  to  e  given  frame  his  bean 
measured  to  be  less  than  1  per  cent  For  the  film  ueeu  to  obteln  the  data,  thle  imojnte  to  *n 
uncertainty  of  *0  8  ft  In  the  object  plan*  The  mtximujj  spetiel  static -resolution  uncertainty 
wet  found  to  bt  >6  ft  The  maximum  dynamic -resolution  uncartiln’y  fella  within  this  iiroti 
The  foreihortening  effect  become#  increasingly  Important  wuh  lncreee#  in  lh#  ehocs-weve 
growth  It  varies  In  *  fixed  manner  ovar  the  range  considered  from  0  1  to  S  5  per  cent  of  the 
radial  dietrnce  measured  Corrections  were  calculated  end  applied  for  each  distance 

Ae  a  result  of  thus  cone  derations  an  average  figure  of  accuracy  of  >3  3  ft  le  assigned  to 
d'.itince  treaejremenle  in  the  tree-air  region 
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3.7  4  Accuracy  of  Preaaure-Dletance  Results 

The  errors  In  Individual  tima  and  distance  measurements  ara  not  carried  over  directly 
Into  tha  praaiura  calculation*  Thai*  arrora  ara  substantially  reduced  by  fitting  the  data  with 
*  curve  by  tha  method  of  lei  at  squares.1  On  the  other  hand,  arrora  In  P,C#  ara  carried  over 
directly  Into  praaiura  calculation*  tt  1*  aatlmatad  that  the  trror  In  that*  data  la  of  tha  order 
of  1  par  cant. 

Tha  error  In  tha  calculated  preeaures  baaed  on  tha  dartvatlvaa  of  tha  fitted  arrival-time 
curvaa  and  tha  atmoapherlc  meeeurementc  mentioned  above  are  conatdarad  to  b«  accurate  to 
1.3  par  cant  at  tha  50-pel  preaaure  level  and  Increasingly  more  accurate  at  tha  higher  levels. 
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CHAPTER  4 

ANALYSIS  AND  DISCUSSION  OF  RESULTS 


In  order  to  compare  the  raritta  obtained  on  King  Shot  -with  those  obtained  on  prevlou#  teat/,, 
It  1*  necessary  to  reduce  the  data  in  tome  coneletent  manner  euch  that  one  variable  le  common 
to  both  aete  of  data.  On  Operation  Turn  jler*  the  Saehe  method*'1  wa#  uaed  to  reduce  the  data  to 
ftandard  eea-level  atmoepherli  condition#,  namely,  14.7  pel  preaeure  and  29S"X  (20*0  tempera¬ 
ture.  In  addition,  the  data  were  i  educed  further  by  applying  the  cube  law  for  the  charge  weight. 
King  Shot  data  have  been  treated  elmllarly.  When  the  data  are  reduced  to  the  equivalent  of  1 
kt(RC)  (radiochemical  kilotonnage)  at  eea-level  condition*,  they  are  termed  “A-ecaled."  Thl* 
term  will  be  uned  frequently  In  the  following  dlecueslon. 

4.1  SCALING  FACTORS  AND  DATA  REDUCTION 

The  acallng  factors  used  to  reduce  the  original  results  contained  L".  Chap.  S  are  given  In 
Table  4.1.  The  arrival-time  data  and  preeeure-dietance  resulte  have  been  A-ecaled.  Aleo,  the 
preesurs-dieiance  results  have  been  reduced  to  sea  level.  The  reduced  result*  have  been  com¬ 
pered  with  these  of  previews  te  le  (1)  to  determine  similarity  of  results,  which  le  actually  a 
test  of  the  scaling  laws,  (3)  to  determine  the  yield  In  terms  of  radiochemical  kllotonnage,  and 
(3)  to  determine  the  TNT  blast  efficiency. 


Table  4.1 —MISCELLANEOUS  DATA  AND  BCALINO  FACTORS 


Assigned  ground  aero  !AGZ) 

N  108,130 

E  124.130 

Actual  ground  aero  (CZ) 

N  108.430  t  30 

X  123.830  «  20 

Burst  height  (h),  ft 

1480  *  20 

Temperature  at  buret  height  (T,),  *C 

23.4 

Radiochemical  yield  (Wrc  ).  kt(RC) 

841  t  30 

Preaeure  at  buret  height  before  shot  (P,).  pel 

13  98 

Factor  to  multiply  preeeuree  to  correct  to  sea  level  (Sp  •  14.7  P,j 

1.033 

Factor  to  multiply  distance  to  correct  to  see  level  [Sc  •  (P,/14.7)s) 

0.9829 

Factor*  to  multiply  distance  to  reduce  to  I  kt(RC)  at  see  level 

0.1206 

[Sj-fP./ll.TWRc^j 

Factor*  to  multiply  time  to  reduce  to  1  kt(RC)  at  eea  level 

0.1217 

18,  -  [IT.  *  *73)  293  Sd) 

A-ecaled  buret  height  ih3d).  ft 

178  *  2.4 

*A-*<al*d  factor# 
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4.2  YIELD  OF  KING  SHOT 
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Before  the  results  ran  be  educed,  a  reasonably  accurate  value  o f  the  yield  (or  radiochemi¬ 
cal  kllotonnage  equivalent  must  be  determined.  A  tentative  value  for  the  yield  was  pjbllshed 
by  Ogle  and  Lolland'  as  550  t  50  kt.  An  Independent  determination  based  on  the  free  air  pres 
sure -distance  results  given  In  Chap.  3  was  made  as  follows 

'.  Values  of  peak  shock  overpressure  were  taken  from  Table  3.3  and  reduced  to  sea  level 
by  using  the  appropriate  scaling  factor. 

2.  These  reduced  pressures  were  located  on  the  A-scaled  pressure-distance  Tumbler 
composite  (Ftg.  4.12  of  reference  1),  and  the  corresponding  Tumbler  distances  were  nut*d. 

3.  The  ratios  of  these  distances  to  the  unreduced  distances  of  King  Shot  given  tn  Table  3.3 
were  determined  and  averaged  over  the  entire  free-air  region. 

4.  The  average  value  obtained  for  the  ratio  of  the  dtstar.cea  was  0. 1 209  t  0  0023.  Equating 
this  value  to  the  A-scaled  factor  (P,  14  ‘IWr^)1'1  and  using  P,  =  IS  98  pst,  the  value  for  Wrc 
was  found  to  be  541  *  30  kt(RC).  This  value  (541  kt)  has  been  used  throughout  the  following 
scaling  procedures 

4.3  COMPARISON  OF  KING  SHOT  DATA  WITH  TUMBLER  COMPOSITES 

4.3.1  Arrival-time  Data 

The  fireball  and  frte-alr  A-scaled  arrival-time  dati  are  given  tn  Table  4  2.  The  free-alr 
data  are  plotted  with  the  A-scaled  Tumbler  composite  arrival-time  curve  in  Fig.  4.1.  The  data 
used  In  formulating  the  Tumbler  composite  curve  and  the  equations  representing  It  are  given 
In  Tables  4.5  to  4.8  and  in  Sec  4  2  1  of  reference  1 . 

The  scale  of  Fig.  4.1,  which  Includes  only  a  small  portion  of  the  Tumblir  composite  curve, 
his  teen  great!"  enlarged  to  nlrn*  the  deviation  of  the  Kl-g  Shot  data  which  fall  within  the  lim¬ 
its  of  msxlmum  uncertainty  of  ths  composits  curvs  The  minimum  uncertainty  In  the  A-scaled 
compos  it*  curve  ts  prlmsrtly  governed  by  the  values  esstgned  to  the  rsdtochemlcal  kHotor-nages 
snd  secondarily  by  the  compound  error*  In  the  ttme  and  spaca  calibrations.  Ths  maximum  un¬ 
certainty  of  all  A-acaled  Tumbler  arrlval-tlms  data  brought  about  by  both  primary  and  second- 
ary  cauaea  wsa  of  the  order  of  2  per  cent.  The  uncertainty  In  the  King  Shot  data  la  of  the  same 
order. 

Only  on  Tumbler  'hots  3  snd  4  were  sufficient  fireball  data  obtained  wltn  which  the  data  of 
King  Shot  could  bs  compared.  1  hia  comparison  la  ahovn  graphically  In  Fig  4  2,  where  good 
agreement  la  tndfated. 

4.3.2  Peax-overpressure  -  Distance  Data 

The  A-acaled  preusure-distance  result*  are  gtven  In  Table  4.3  snd  are  shown  grsphlcal’v 
In  Fig.  4.3  The  A-scaled  Tumbler  composite  pressure  distance  curve  1*  tncludel  for  com¬ 
parison.  The  data  used  in  the  formulation  of  the  Tumblsr  composite  curvs  can  be  ‘ound  tn 
Tables  4.11  lo  4.14  of  reference  l 

On  ths  sversge,  the  A-scsled  King  p. enure -dietin  '#  results  agree  with  the  Tumbler  com¬ 
posite  to  within  2  0  per  cent.  At  the  extremes  of  the  pressure  rsr.gr  measured,  however,  the 
Kmg  dsts  are  low  by  approximately  8  per  cent  Figure  4  3  represents  the  best  possible  over¬ 
all  fit  with  ths  Tumbler  composite  snd  provtdei  the  basis  for  the  yteld  calculation  of  Sec.  4.2 

4.4  TNT  BLAST  EFFICIENCY 

Bv  fitting  TNT  pressure  distance  data  to  dal*  from  King  Shot,  s  vslue  for  the  TNT  effi¬ 
ciency  of  the  nuclear  weap  mi  can  be  delerm.ned  The  procedure  used  Is  is  fellow* 

'  The  free  air  preasure  distance  data  for  the  nuclear  expl'-slcn  found  in  Table  3  3  were 
c  rrurted  to  aea  level  Table  4  3  and  are  plotted  In  Fig  4.4, 
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Tabl*  4.2— A-SCALED  ARRIVAL-TIME  DATA 


Distance  from 
weapon  zero, 
A-scaled,  ft 

Arrival  time, 
A-scaled, 

sec 

Distance  from 
weapon  zero, 
A-scaled,  ft 

Arrival  time, 
A-sceled, 

sec 

70.1 

0.00133 

288.1 

0.04429 

91.8 

0.00261 

291.7 

0.04572 

108.3 

0.00424 

294.6 

0.04715 

120.5 

0.00567 

299.8 

0.04858 

132.0 

0.00710 

304.3 

0.05001 

141.8 

0.00853 

307.5 

0.05144 

180.9 

0.00996 

311.3 

0.0  5  287 

159.1 

0.01139 

314.9 

0.05430 

a  a  a 

e  e  e 

318.8 

0.05573 

177.2 

0.01425 

323.4 

0.05716 

184.3 

0.01568 

325.5 

0.05859 

190.4 

0.01711 

329.6 

0.06002 

197.0 

0.01854 

331.6 

0.06145 

202.7 

0.01997 

336.6 

0.06  288 

210.3 

0.02140 

338.2 

0.06431 

217.0 

0.02283 

343.2 

0.06574 

221.2 

0.02426 

345.7 

0.06717 

226.1 

0.02570 

349.7 

0.06860 

231.5 

0.02713 

352.3 

0.07003 

«*•*  • 

89V,« 

0.0285* 

354.8 

0.07146 

244.0 

0.02999 

358.8 

0.07289 

24C.8 

0.03142 

362.2 

0.07433 

252.0 

0.03285 

366.3 

0.07576 

257.8 

0.03428 

368.8 

0.07719 

263.3 

0.03571 

373.0 

0.07861 

265.8 

0.03714 

375.3 

0.08005 

270.7 

0.03857 

379.3 

0  08148 

273.5 

0.04000 

381.3 

0.08291 

279.3 

0.04143 

.  a  . 

•  •  e 

282.9 

0.04286 

387.4 

0.08578 

390.8 

0.08721 

2.  The  theoretical  A-scaled  TNT  pressure -distance  data  presented  in  Table  4  4  end  shown 
In  Fig.  4.5  are  baaed  on  Information  published  by  Kirkwood  and  Brinkley'  and  Hermann  *  For 
equivalent  preeauraa  the  corresponding  dittance*  were  read  from  th»  curves  of  Fig.-,  4  4  end 
4.5.  The  ratio  of  the  nuclear  distance  to  the  TNT  dletance  for  each  pressure  level  is  equal  to 
the  cube  root  of  the  TNT  kllotonnage  [ktITNT!)  equivalent  et  that  preeejre. 

3.  The  TNT  efficiency  for  each  preeeure  level  was  determined  by  cubing  the  distance  ratio 
and  dividing  by  the  radiochemical  yield  of  the  nuclear  weapon.  The  average  efficiency  via  taken 
over  e  given  preee  re  range,  the  upper  1'mlt  of  which  le  of  the  order  of  200  pel.  because  f  r 
preeeuree  higher  than  th.e  '.he  elope  of  the  TNT  preeeure-dietance  curve  falls  off  rapidly  e§ 
compared  to  that  for  a  nuc.ear  explosion  (see  Fig.  4.5). 
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Table  4.2 —  SCALE’.)  PEAK-OVERPRESSURE -DISTANCE  DATA 


Distance  from 

Distance  from  weapon 

Distance  from 

Peak  overpressure, 

weipon  *<ro, 

ssro,  reduced 

weapon  ssro, 

reduced  to  set  level, 

ft 

to  sts  level,  ft 

A-scaled,  ft 

pel 

000 

390 

72 

3086 

700 

MS 

14 

3116 

800 

788 

68 

3034 

900 

885 

101 

1454 

100u 

983 

111 

1070 

1100 

1081 

113 

623 

1200 

1179 

146 

845 

1300 

1278 

161 

622 

1400 

1276 

166 

438 

1300 

1474 

161 

136 

1600 

1673 

162 

268 

1700 

1871 

206 

231 

1800 

1789 

*17 

214 

«90C 

1888 

226 

181 

3000 

1988 

241 

136 

2100 

2084 

263 

136 

3200 

2182 

265 

132 

3300 

2261 

277 

107 

2400 

2339 

289 

67.0 

2500 

«*af  1 

301 

67.0 

2000 

2338 

314 

78.8 

2700 

2634 

336 

70.9 

2600 

2732 

331 

63.2 

2900 

2830 

350 

57.7 

3000 

2949 

362 

52.3 

3100 

3047 

374 

47.7 

3200 

2143 

386 

43.2 

4.  The  average  value  (or  the  TNT  efficiency,  whe..  mul'  piled  by  the  radiochemical  yield  of 
the  nuclear  weapon,  glvee  a  value  for  the  yieid  In  terms  of  TNT  kUotonnage.  Uee  of  this  value 
In  the  King  scaling  operation  produced  the  results  shcrvn  In  Elg.  4.3  In  which  the  fit  to  the  TNT 
data  Is  Illustrated. 

Following  this  outline,  the  average  TNT  efficiency  of  King  Shot  was  found  to  be  38.8  t  2.0 
per  cent  which  corresponds  to  a  yield  of  309.0  Irt(TNT)  over  the  pressure  range  cf  200  to  30  pel. 

The  percentage  TNT  efficiency  as  obtained  above  Is  lows.-  than  the  Tumbler  composite 
average  of  42  1  per  cent  within  the  200-  to  30 -pel  pressure  range.  This  Is  due  chiefly  to  the  fall- 
off  of  pressure  at  the  lower  end  of  the  pressure  range  measured.  It  should  be  etreseed  again1  * 
that  the  value  of  TNT  efficiency  Is  vary  aervltlYS  lo  small  differences  In  the  distance#  at  whir*- 
e^ul'al;  pressures  sre  determined  frem  one  shot  to  the  next.  Becaute  cf  the  cube  lav  th.t  Is 
brought  Into  the  cs’culatlons.  these  small  differences  tre  significantly  magnified.  Caution  must 
be  exercised  when  Judging  different  sets  of  dais  on  this  basis. 
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Table  4.4  —  COMPARISON  OF  KING  SHOT  DATA  AND  TNT  DATA 


Peak  over¬ 
pressure, 
pal 

Radius  for 

1  kt  of  TNT 
at  sea  level, 
ft 

Distance  from 
weapon  aero, 
reduced  to  sea 
level,  ft 

TNT 

equivalent 
of  King, 
kt(TNT) 

TNT  efficiency 
of  King,  % 
[kt(TNT)  x  100] 

[  kt(KC)  ] 

Distance  from 

weapon  aero,  re¬ 
duced  to  1  kt(TNT) 
at  sea  level,  ft* 

eoo 

132 

184 

700 

145 

194 

600 

160 

205 

500 

179 

218 

400 

204 

237 

300 

241 

264 

200 

295 

1815 

232.8 

43.0 

306 

150 

334 

2015 

219.5 

40.6 

339 

100 

393 

2340 

211.1 

39.0 

394 

90 

410 

2435 

209.5 

38.7 

410 

80 

430 

2540 

208.0 

36.1 

428 

70 

413 

2670 

202.0 

37.3 

449 

60 

480 

2810 

200.6 

37.1 

473 

50 

SIS 

299  3 

196.7 

36.4 

504 

Av. 

209.6  a  10.7 

38.8  a  2.0 

*Valuaa  scaled  on  ths  basis  of  209.8  kt(TNT). 
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CHAPTER  5 

CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  INSTRUMENTATION 

The  photographic  results  wcr*  good  daaplta  tha  facta  that  ont  Camara  tailed  to  operate, 
that  poor  atmoapherlc  condltlona  azlated,  and  that  timing  marks  were  applied  at  a  low  rate. 
Timing  algnal  ralaya  and  a-c  powar,  houaad  In  a  new  and  vary  much  chaaper  aheltar,  worked 
properly,  and  all  IS  rcckete  fired  aa  planned.  The  fan-type  grid,  ueed  because  of  tha  extreme 
apace  Umltatlona  Imposed,  proved  to  be  very  satisfactory.  Compared  to  tha  vertical-line  grid 
ueed  on  previous  testa,  three  disadvantages  were  noted. 

1.  Mors  rockets  must  be  fired  In  order  to  achieve  equivalent  space  coverage. 

I.  The  rocket  battery  must  be  reasonably  compact,  which  makea  the  loadtng  operation 
more  hazardous. 

3.  The  probability  of  a  trail  coinciding  with  a  radial  line  from  the  burst  Is  Increased.  This 
la  undesirable  because  the  light  refracted  from  such  a  trail  cannot  be  readily  detected. 

As  a  result,  It  la  recommended  that  the  fan  grid  be  used  only  If  necessary. 

5.3  DATA  ANALYSIS 

The  use  of  IBM  equipment  decreased  the  time  required  to  fit  an  analytical  expression  to 
the  measured  arrival-time  data  by  35  per  cent.  Identical  results  were  obtained  using  both  the 
new  and  the  old  methods;  thus  complete  confidence  can  be  placed  In  all  the  results,  old  and  new 
alike.  The  average  deviation  of  the  data  points  from  the  fitted  curve  was  lees  than  1  per  cent 
In  distance  or  time. 

5  3  RESULTS 

Arrival-time  data  were  obtained  throughout  the  first  0.7  sec,  during  which  time  the  shock 
wave  traveled  out  to  approximately  3350  ft  Baaed  on  these  data  the  following  results  were  ob¬ 
tained  for  King  Shot- 

I.  The  yield  of  54!  t  30  kt(RC)  determined  from  the  comparison  of  King  results  with 
Tumbler  composite  free-alr  pressure  reaulte  te  In  good  agreement  with  the  value  of  550  »  50 
ktiRC)  given  by  Ogle  and  l  aland.1  On  this  bails  the  range  over  which  the  seeling  lews  for  free 
elr  can  be  ueed  is  extended  to  550  kt(RC). 

3.  The  TNT  blast  efficiency  of  King  was  of  the  earns  order  as  that  of  the  Tumbler  average. 
38  8  and  43.'  per  cent,  re  tpectlvely,  in  the  preeaure  range  of  300  to  50  pat.  It  ta  interesting  to 
note  that,  aa  the  kilotonnage  Increases,  the  TNT  efficiency  apparently  decreases  slightly  > see 
reference  3  for  comparleon). 
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S.  In  the  fireball  region 

R  ■*  3302. 3t*  ***  R  s  900  ft  (5.1) 

P  a  R“,  w 

whereas  In  the  free -air  region 

t  *  o.ooMi  -  y  ioslaiT8;1"!-.  anj  -  0.6236  900  ft  *  r  s  3350  ft  (3.2) 

where  R  la  the  radial  dletance  from  weapon  zero  In  feet,  t  la  the  tiros  In  seconds,  and  P  is  ths 
peak  shoca  overpressure  In  pounds  per  square  Inch. 

5.4  THERMAL  EFFECTS 

The  precursor  wave  detected  by  the  pressure-time  gauges  of  Project  3.1  ovsr  ths  island 
at  Run'.t  could  not  be  observed  In  the  Project  6.13  films  because  the  view  at  the  Island  was  ob¬ 
structed  by  palm  trees  on  an  Island  in  the  foreground,  and  the  camera  location  was  poor  in  this 
respect.  This  situation  was  not  remedied  before  the  teat  because  of  insufficient  tiros  and  man 
power.  In  the  event  of  a  similar  taat  In  the  f  iture,  every  effort  should  be  mads  to  assure  an 
unobstructed  view  of  GZ. 

The  films  did  not  show  a  precursor  wave  over  ths  water,  although  the  view  was  unob. 
structed.  This  confirms  the  records  produced  by  the  gauges  placed  ovsr  ths  water  by  Project 
6.1.  These  records  showed  that  the  shock  wave  had  a  fast  rise  tims  and  decayed  In  ths  manner 
expected  (see  Ivy  Project  6.1  report,  WT-603) 

No  additional  conclusive  Information  was  obtained  to  explain  the  deviation  of  Gresnhouss 
free-air  pressure-distance  results  from  the  Tumbler  composite.  There  is  only  one  thuught 
concerning  this  problem  wherein  the  King  Shot  results  might  Indicate  a  solution.  On  the  basis 
at  relative  energy  Incident  at  the  ground,  Greenhouse  Easy  was  more  than  twice  as  effectivs  as 
King;  yet  a  precursor  wave,  believed  to  be  formed  as  the  result  at  the  existence  at  a  thermal 
layer, ,_4  was  observed  on  King  Shot.  It  is  thus  possible  to  conceive  that  ths  atmospheric  heat¬ 
ing  associated  with  a  thermal  layer  on  Greenhouse  was  so  Intense  that  the  shock  wave  In  fres 
air  traveled  much  faster  than  it  would  In  an  unheated  atmosphers.  Since  temperature  meas¬ 
urements  observed  Just  prior  to  the  burst  on  Greenhojss  were  used  to  compute  'he  some  ve¬ 
locity  (which  Is  used  In  the  pressure  calculations,  Eq  1.5),  It  is  also  conceivable  that  the 
values  used  were  significantly  low.  This  would  explain  ths  high  pressures  obtained  at  Green¬ 
house  as  compared  with  those  from  Tumbler,  at  least  in  part.  The  proximity  of  the  ground  and 
other  arguments  proposed  In  references  3  to  4  may  also  have  contributed. 
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APPENDIX 

THE  METHOD  OF  REDUCING  THE  DATA 

By  T.  8.  Walton 


A.l  NATURE  OF  THE  DATA 

The  experimental  data  conalat  tn  a  larga  number  of  pair*  of  poaltlon  and  ttma  maaaura- 
manta  obtained  from  a  high-speed  motion -picture  film  of  tha  shock  wave  from  an  a/.ploalon. 
The  problem  involved  In  reducing  the  data  le  baalcally  that  of  finding  a  mathematical  function 
which  effectively  correlatee  theae  obaerved  valuea  of  poaltlon  end  time,  ao  that  the  velocity  of 
the  ahock  front  can  be  determined  from  the  derivative  of  tho  function. 

Some  method  of  fitting  the  data  which  le  baaed  on  the  principle  of  lee  it  aquiree  aeeme 
appropriate  becauae  the  meeaureminta  are  not  exact  but  contain  preeumably  random  errora 
from  many  eourcaa.  The  ehock  wave  ltaelf  may  not  be  propagated  with  perfect  regularity 
becauae  of  the  alight  lnhomogenelty  of  the  atmoaphere,  and  aalda  from  thla  there  are  numeroua 
errora  introduced  by  the  meaeurtng  equipment  (for  example,  fluctuation!  In  the  tyeed  of  the 
camera  motor  and  uncertainty  in  the  location  of  the  ahock  front  due  to  the  optical  resolution 
of  the  lane  eyatem  or  the  grain  elte  of  the  film). 

Prior  to  the  completion  of  the  analyele,  It  could  not  be  ascertained  whether  the  dlsperelon 
of  the  data  was  primarily  the  reeult  of  errors  In  the  position  or  the  time;  so  attempts  were 
made  to  fit  the  observations  both  ways.  In  order  to  keep  the  snalysls  simple  and  minimize  the 
amount  of  computation,  polynomials  In  the  time  and  In  the  distance  were  first  tried  and  then 
later  the  ratio  of  two  simple  polynomials.  However,  these  were  generally  unsatisfactory  be¬ 
cause  they  would  fit  only  limited  stretches  of  dsta.  The  Junction  points  between  succeeelve 
curves  were  not  smooth  and  showed  very  abrupt  changes  In  elope. 

The  entire  range  of  observations  couia  have  been  approximated  io  any  desired  degree  by 
Introducing  many  more  arbitrary  patsmeters;  but  In  such  a  csss  ths  approximating  function 
would  have  too  much  "flexibility,''  and  Its  derivative  could  oscillate  wildly  throughout  the  range 
of  data  points,  leading  to  meaningless  results. 


A.J  THEORETICAL  BA3U  FOR  THE  ANALYSIS 

The  limited  success  achieved  with  arbitrary  polynomials  indicated  that  a  theoretical  basis 
for  selecting  an  approximating  function  should  be  sought.  It  seemed  desirable  to  find  i  form  of 
fundi  -n  which  corresponds  a«  nearly  ea  possible  with  the  actual  physical  behavior  of  an  ex 
plosion  and  also  one  In  which  the  arbitrary  parameters  represent  degrees  of  freedom  having 
phyelca,  counterparts  The  reliability  of  the  fit  could  then  be  checked  not  only  by  observing  the 
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trend  of  the  residuals  but  also  by  comparing  tha  a  poatarlorl  values  of  tha  parametera  with 
reasonable  a  priori  estimates  of  them. 

A  simplified  approach  to  tha  problem  la  to  treat  tha  air  surrounding  the  blast  as  a  homo¬ 
geneous  medium  of  Infinite  extent  obeying  the  aquation  of  state  of  a  perfect  gas.  To  Idea. ltze 
the  proolem  further,  It  is  assumed  that  all  the  energy  la  suddenly  released  at  a  point  within 
the  medi-iti  (as  heat  energy  available  for  performing  work  on  the  air  but  none  being  radiated!. 
Then  an  attempt  is  made  to  deduce  from  the  theoretical  equations  of  gas  dynamics  the  law 
which  describee  the  motion  of  the  shock  wavs  and  tha  manner  In  which  It  depend?  on  the  five 
parameters  Involved;  namely,  the  original  density,  tsmperatura,  and  specific  kea'  ratio  of 
tha  air  and  tha  total  energy  and  time  of  Inception  of  the  blast. 

A  ssarch  of  tha  lltsraturs  revealed  several  noteworthy  articles  dealing  with  this  problem. 
The  basic  theory  wee  first  developed  In  1941  by  Taylor,1  the  end  result  being 


t  -  c 


1 

U 


£1 

dR 


(A.l) 
(A.  2) 


where  R  •  radii  a  of  the  shock  front 
U  ■  velocity  of  tha  shock  front 
t  -  time  of  arr‘  -el  of  the  shock  at  R 
c  ■  time  at  which  the  blast  began 
C  -  total  anargy  Instantaneously  released 
K  -  s  parameter  depending  on  the  specific  best  rstlo 
p  ■  original  density  of  the  air 

The  only  additional  assumption  implicit  In  Taylor's  solu.lon  Is  the  strong  shock  condition, 

I  s.,  thit  the  Rankins -Hugonlot  relations  assume  their  asymptotic  forms  corresponding  to  an 
Infinite  pressurs  ratio  across  the  shock  front.  TMs  assumption  Is  correct  ss  long  ss  the  Mach 
number  of  tha  shock  la  quite  Urge,  but  It  is  clearly  wrong  In  the  limit  as  R  Increases  indefi¬ 
nitely;  for  according  to  Eq.  A.2  the  velocity  of  the  shock  will  ultimately  drop  to  xero  (rather 
than  approach  tha  velocity  of  sound,  as  It  must).  Taylor’s  formula  Is  correct  for  large  radii 
only  when  the  temperature  of  the  medium  Into  which  the  shock  sxpands  Is  at  absolute  zero,  so 
that  both  the  preseure  and  tha  vsioclty  of  sound  can  ba  aaro  outside  while  the  density  remains 
finite. 

Soma  progress  toward  removing  this  shortcoming  In  Taylor’s  solution  has  recently  been 
made  by  Newton,1  who  modlflse  toms  of  tha  flow  aaeumpttons  to  allow  for  tha  raqulred  bs- 
havlor  al  both  large  and  small  raan.  Nswton  shows  that  his  equations  do  not  violate  try  physi¬ 
cal  principles  but  he  Is  unabls  to  obtain  a  completa  solution  bacauaa  of  tha  great  complexity 
of  tha  equations.  However,  he  doea  find  limiting  aoluilons  for  both  small  and  large  radii.  The 
reault  aa  R  approaches  aaro  'grass  with  that  given  by  Taylor,  whils  for  larga  rad!)  he  finds 
that  the  velocity  of  the  shock  wava  exceeds  that  of  a  sound  wavs  by  an  amount  which  falls  off 
aa  1/R1,1.* 


*  This  decay  law  for  the  excess  velocity  of  spherical  shock  waves  la  somewhat  at  variance 
with  the  result*  of  Brinkley  and  Kirkwood1  and  Whilham.*  They  obtain  the  formula  1  R  [.  'g 
(R  P-elj*  V  where  R,  It  some  suitable  constant.  However,  In  each  cea*  their  concluaion*  are 
arrived  at  on  the  assumption  that  at  autflcler.tly  great  distances  the  shock  wevt  becomes  sub- 
eta  lally  an  acoustic  wav*  end  thal  the  flow  acroea  the  ahock  front  can  then  be  considered 
leentropic  By  contraet,  Newton  use*  the  sxnct  Ranklns-Hugonlot  relatione,  end  the  change  in 
r  i  op>  although  seemingly  negligible,  maybe  the  crux  of  the  matter. 
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Newton’a  reeult  immediately  auggeete  a  eimple  type  of  algebraic  function  which  might  be 
ueed  to  appro:  mat*  tha  valoclty  of  tha  ahock  for  all  radii,  namely, 


l\ 

rt 

s) 

(A.S) 

where  a  la  the  velocity  of  an  acouatlc  wave  In  the  region  ahead  of  the  ahock  and  b  la  a  ecale 
paramater  which  dependa  on  the  Initial  denalty  and  temperature  of  the  air  and  the  amount  of 
energy  releaaed.  U  thla  relation  la  aeaumed  to  be 


b* 


then  Eq.  A.S  can  be  written 


U  -  a  ♦ 


RXr 


(A.  4) 


(A.  5) 


and,  If  a  la  neglected,  thla  reducea  to  Eq.  A.S.  Thua  Taylor'a  formula  for  the  velocity  haa 
merely  been  augmented  by  the  amount  a  to  account  for  the  finite  -eloclty  of  propagation  of  the 
wave  aa  R  approachea  Infinity,  and  thla  la  compatible  with  hia  aolutlon  aa  R  approachea  aero 
alnca  In  that  caae  the  aecond  term  on  the  right  aide  of  Eq.  A.S  becomea  Indefinitely  large  and 
the  added  contribution  of  the  fir  at  term  la  then  lnatgnlftcant. 

If  reaaonable  agreement  with  a  aet  of  experimental  data  can  be  obtained  by  adjuatlng  the 
parametere  a  and  b  In  Eq.  A.S,  thla  quaal -theoretical  formula  can  be  aubatrntlated.  In  order 
to  carry  out  an  actual  fit,  It  la  neceaaary  to  determine  the  mathematical  relation  between  time 
and  dletance  which  thla  formula  lmpllea.  The  time  of  arrival  of  the  ahock  at  the  point  R  la 


«~jr 

Subatltutlng  the  expreaalon  (Eq.  A.S)  for  U, 


(A-fl) 


(A.T) 


la  obtained.  For  convenience  in  evaluating  the  Integral,  Eq.  A.7  la  written  aa  followe: 


t-  c  ♦ 


where  a  •  (Ttb )*'*.  Equation  A.S  may  ba  Integrated  to  give 


Hf  •(*) 


de 


(A.S) 


R  b  L  n  a/T  a  ♦  1  1 

c  ♦ - —  S/S  arctan  - - log  -f- - L- 

a  Sa  l  J-a  a*  -a  ♦  1 


(A.S) 


The  aum  of  the  fir  at  two  terma  on  tha  right  aide  of  Eq.  A.S  correepcnde  to  the  time  of 
arrival  of  an  acouetic  wave  at  the  position  ft,  ard  the  iaet  term  eccour.te  for  the  lead  time  of 
t  epherlcel  ehock  wave  ever  an  acouetic  wave.  It  te  InUreetirvg  to  note  that  the  tn'.egral  oc¬ 
curring  In  Eq  AS  haa  an  aeymptotlc  value  of  2.4184  when  »he  u^per  limit  of  Integration 


4$ 


approach**  infinity.  Thia  can  ba  lntarpratad  by  saying  that  a  spherical  ehock  wava  wilt  over¬ 
take  ar  acouattc  wave  which  previously  emanated  from  tha  earn*  point  in  apaca,  provided  that 
tha  latter  did  not  have  a  head  atart  exceeding  2.4184  fb  a)  tima  unit*.* 

Each  of  the  parameter*  in  Eq.  A. 9  ha*  a  phyalcal  Interpretation,  and  thia  ahould  facilitate 
the  compariaon  of  fitted  data  obtained  from  different  teat*.  The  quantity  c  la  aimply  the  time 
intercept,  a  la  the  asymptotic  value  of  the  velocity,  and  b  la  proportional  to  tha  phyalcal  di¬ 
mensions  of  the  blast,  which  vary  as  the  cubs  root  of  the  energy  released.  The  well-known 
scaling  laws  for  explosive  phenomena  may  be  applied  to  Eqs.  A.S  and  A. 9.  For  example,  the 
Mach  number  of  the  shock  fU/a)  depend*  only  on  the  nondlmsnelonal  r  .tto  (R/b),  and  con¬ 
sequently  the  trajectory  of  the  shock  front  will  be  similar  for  an  exp'  >*lon  of  any  size  in  a 
given  medium  If  the  time  and  distance  ecalaa  are  both  multiplied  by  .  factor  proportional  to  b. 

Thu*  It  is  seen  that  the  data  may  be  praaanted  In  terme  of  any  da  tired  units  of  distance 
and  time.  The  time  values  may  also  contain  an  additive  constant,  1  e.,  zero  clock  time  need 
not  correspond  to  the  start  of  the  blast.  However,  the  radial  dtstancse  must  be  referred  to 
the  true  center  of  the  explosion,  although  they  may  be  ecaled  by  an  arbitrary  factor. 


A.S  COMPUTATIONAL  PROCEDURE 


In  this  section  Is  given  an  outline  of  a  procedure  for  determining  the  valute  of  the  param¬ 
eters  'n  the  approximating  function  by  the  method  of  least  squares.  Although  In  principle  the 
process  of  fitting  a  given  set  of  data  could  be  carried  out  by  manual  computation,  It  would  be 
very  laborious  and  time  consuming;  so  It  will  be  presumed  that  the  routines  described  In  the 
following  are  to  be  programed  for  automatic  digital  computing  machine*. 

Let  the  symbol  f  «  f[a,b,c,R]  stand  for  the  mathematical  function  which  approximate* 
the  measured  time  tjft],  the  quantities  In  bracket*  Indicating  the  functional  relations  Involved. 
For  convenience,  Eq.  A. 7  la  rewritten  a* 


where 


t  -  c  «  C*  -  T)  - 
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Thia  permit*  the  calculation  to  be  carried  out  In  terme  of  nondlmtnslonal  quajitlc:  •  »  x  and  y. 
To  almpllfy  the  analysis,  1!  Is  aseumed  that  me  value*  of  R  ar*  given  exactly  ard  that  all  the 
value*  of  t  hava  been  corrupted  by  (tails. leal  arrera  of  uniform  dlaperalon.  The  validity  of 
thia  aaaumptlon  can  be  checked  by  a  caraful  inspection  of  the  final  result*  of  the  computation 
The  reatdual  error  associated  with  the  1th  data  point  la  daflned  aa 


* 


i 


-t. 


(x1-y1)-*c-t, 


(All) 


*  Bv  way  of  contrast,  tha  decay  law  discussed  In  th*  preceding  footnote  give*  a  lead  time 
proportional  to  th#  fsetor  •'Tog  iR  R,)  which  has  no  limit  ts  Ft  approaches  infinity  This  In- 
<“. rates  that  a  spherical  shock  wav*  (regardless  of  how  weak)  would  evantuaJly  overtake  any 
acoustic  wave  regard. taa  of  how  much  eerlier  11  haa  been  emitted,  a  result  contradictory  to 

Cuti.n.of,  ar  se 
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And  to  make  the  problem  definite  U  1*  auppoaed  that,  In  accordance  with  the  theory  of  proba¬ 
bility,  the  "beet"  eet  of  vaiuee  for  the  parameter!  a,  b,  and  c  hat  been  found  when  the  eurr.  of 
equaree  of  ail  the  reeiduale  la  a  minimum.  Now  the  reelduai  equation  !Eq.  A.  11)  la  seen  to  In¬ 
volve  the  parameters  1/a  and  c  In  a  linear  manner ,  but  thla  le  not  true  of  b  on  account  of  the 
wey  in  which  a  and  y  depend  on  b.  Consequently  the  uaual  procedure  for  obtaining  the  normal 
equatione  for  determining  the  beet  valuea  of  the  parametera  would  lead  to  an  Intractable  eye- 
tern  of  nonlinear  equatione  (with  no  unique  aolution  but  an  Indefinite  number  of  aolutiona,  among 
which  not  more  than  one  wou'd  be  aubject  to  phyatcal  Interpretation). 

To  avoid  thla  difficulty,  tha  method  of  “differential  correction#"'  for  iteratively  improving 
an  Initially  (uaaaad-at  aat  of  valuea  la  uaad.  Tha  function  f  defined  by  Eq.  A.  10  can  be  expanded 
In  e  Taylor  aerlea  about  tha  point  (a,b,c,Ri).  Thua  tha  change  induced  In  f  by  email  changea  in 
a,  b,  and  c  may  be  expressed  aa 

Of  «  da  ♦  ob  ♦  Ac  ♦  higher  order  term#  (A.  12) 


Provided  that  tha  initial  valuea  are  choaen  efficiently  doat  to  the  deelred  eolutlon,  all  terme 
beyond  the  flrat  order  may  be  neglected,  and  the  reelduai  aatoclaied  with  the  fth  point  becomee 


(fj  ♦  Aft)  -  t,  •  (f,  - 1 p  ♦  A ,0a  ♦  B,Ab  ♦  CtOc  (A.13) 

where  At  •  (a^/aa)  hold!.  g  b,  c,  and  Rt  conatant;  B,  «  (9ft/0b)  holding  a,  c,  and  R,  conatant;  and 
C(  ■  <»f,/»c)  holding  a,  b,  and  Rt  conatant. 

Tha  problam  now  reduce!  to  finding  tha  correctlone  which  make  tha  eum  o I  equarea  of  tha 
new  raalduala  a  minimum.  Thla  lead!  to  a  eyeteca  of  three  aimultaneoue  linear  aquatlona  In  tha 
unknown!  Aa,  Ab,  and  Ac,  namely, 

(EA*)Aa  ♦  (TAjB^Ab  ♦  (IAjC^AC  -  lA^t,-!,) 

(IBiAilAa  ♦  (IBj)Ob  ♦  (IB^JAc  -  IB^-f,)  (A.14) 

irC(A,)Aa  ♦  (ICE^Ab  ♦  (EC*) Ac  -  EC^-f,) 


Tha  aolution  of  tha  ayatem  (Eq.  A.14)  glvaa  tha  ftret-order  correction!  to  be  added  to  the  Initial 
valuea  of  a.  b,  and  c.  Tha  reeulia  may  be  Improved  by  repeeting  the  proceee,  uelng  the  "cor¬ 
rected"  vaiuee  In  place  of  the  original  choicee.  l'owever,  thte  method  will  not  alwey*  converge 
to  the  deelred  eolutlon  if  the  initial  gueeeee  ar*  r,  •  accurate  enough  or  if  the  dleperelon  of  the 
data  due  to  random  error*  la  aaceaalva. 

An  appropriate  initial  value  for  a  1*  to  ua*  tha  vtloclty  of  ecund  corrvapondtng  to  tha  p.-e 
vailing  etmoaphertc  condition#  at  tha  time  of  tha  blaat.  Tha  value  of  b  might  be  estimated  on  the 
bail*  of  the  equivalent  energy  of  the  bleu  (If  thte  were  known)  a*  indicated  in  Eq.  A. 4.  Alterna¬ 
tively  the  quantity  0  IS  **0*  can  be  subetituted  for  Kp/E  In  Taylor'*  eolutlon  for  email  R  (Eq. 
A.t(,  the  result  being 


(t-cj* 


0  1SR* 


( A.  IS) 


Thte  can  also  be  obtained  by  dropping  the  eddltlv*  conatant  1  In  Eq.  A. 3  beiore  carrying  out  the 
integral. on  of  1  U.  For  the  ith  date  point,  Eq.  A.  13  may  be  written 

b  * 


(A  IS) 
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Equation  A.  1C  can  bt  jaed  to  aatlmata  b  provided  that  a  pair  of  data  vaiuaa  R(  and  t(  art  taken 
sufficiently  near  the  origin  of  the  bit  it  that  the  Mach  number  of  the  ahock  at  that  point  it  targe 
compared  to  unity.  Obvloutly  thli  procedure  require*  eattmatej  of  both  a  and  c,  whereat  the 
procedure  bated  on  the  equivalent  energy  of  the  blaat  doet  not  require  any  knowledge  of  c.  How¬ 
ever.  the  origin  of  the  time  icale  can  ueuaily  be  found  from  Independent  mtaiurementi. 

Aalde  from  the  foregoing,  no  preliminary  value  la  needed  tr>  the  caae  of  c  alnce  It  entera 
the  definition  of  f  (Eq.  A.  10)  in  >  completely  linear  manner,  and  it  followa  that  the  solution  of 
the  ayttem  of  Eq.  A.  14  la  independent  of  c.  Thla  la  readily  aeen  by  examining  the  expretalont 
uted  for  calculating  the  valuea  of  At,  B(,  and  C(.  Thua,  carrying  out  the  partial  dlfftrentlatlona 
of  f  with  reapect  to  a,  b,  and  c  aa  Indicated  following  Eq.  A.lS, 

A1  ”  ~  (^)  “  Tl> 

B‘ '  (1)  (rrlrt  -  t,)  (A.n, 

c,-l 


The  required  valuea  of  the  definite  Integral  y  can,  of  couree,  be  obtained  from  tablet  of  the 
lnverie  tangent  and  natural  logarithmic  function!  aa  Indicated  in  Eq.  A.9.  However,  when  the 
computation!  are  to  be  carried  out  on  an  automatic  digital  computing  machine,  the  entire  pro¬ 
cedure  can  be  greatly  expedited  by  evaluating  y  from  an  analytic  continued  fraction,  namely, 


-i 


1  ♦! 


T1 


■J.  r 


♦  t* 


,  4t« 

9a> 


t  ♦ 


11  ♦ 


S9a> 


14  ♦ 


04a1 


17  ♦ 


lit 
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(A.18) 


where  a*  -  x  ■  R  ’b.  Thie  it  a  epeclal  cate  tfor  the  1.5-power  index)  of  a  formula  given  by  Wall.* 
A  doien  or  to  termt  of  thla  continued  fraction  art  tuffident  to  determine  accurate  valuea  of  the 
Integral  for  argument!  up  to  t  -  2  (1.*.,  f<  r  radii  not  exceeding  4b),  at  which  point  the  velocity 
of  the  ehock  has  dropped  to  1.125  time*  a.  Although  it  it  not  iiiriy  to  occur  In  practice,  • 
greater  range  for  the  upper  limit  could  be  handled  by  employing  more  terme  alnce  the  continued 
fraction  (Eq  A.  IS)  converge*  for  all  poiltlve  valuea  of  l. 


A. 4  CONCLUSION 

The  reeulte  obtained  by  fitting  aevaral  eete  of  data  to  a  formula  of  the  type  deecrlbed  inJi- 
cete  that  Eq  A.9  la  a  gxid  approximation  to  tha  trajectory  of  the  ehock  from  e  etrong  expl  ’»ion 
The  reeiduai*  for  eome  of  the  moet  tellable  data  ahow  a  d.eper*'on  of  about  2  maec  In  the  time 
reran  .ire  men!*  over  the  entire  range  of  daia.  However,  they  are  noi  completely  rarJim  but  ex¬ 
hibit  a  number  of  p<  altlve  and  negative  grouping*  The  pnnod  of  this  fluciuiiion  dure  .not  appear 
U  charge  »,th  increaaing  radiu*  but  remain*  practically  cor.etani,  which  euggeat*  that  the  acil- 
1  at  ton  II  nut  Inherent  in  the  blaet  phenomenon  but  results  Inatead  from  a  deficiency  In  the  meuh 
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od  of  meaaurlng  the  time  (for  tumplt,  variations  In  camera  speed).  There  le  alto  tome  evi¬ 
dence  that,  if  thlt  periodic  oecllLation  couid  be  taken  out  of  the  time  data,  the  retnaii  lug  diaper- 
tlon  could  be  attributed  to  random  errort  of  about  2  ft  In  the  location  of  the  ehock  front  deter¬ 
mined  from  each  frame  of  film. 

From  a  ttudy  of  the  data  available  at  thlt  tlmt,  It  tetmt  beat  to  ascribe  the  Inherent  er¬ 
rors  to  the  tlmt  measuremente,  and  thlt  leade  to  the  simplest  method  of  treatment.  If  It  should 
ever  prove  desirable  to  ateume  that  the  timet  are  given  accurately  arid  that  errort  are  present 
In  the  positions  only  with  a  urlform  dispersion,  It  would  (till  be  possible  to  employ  the  proce¬ 
dure  described  here.  All  that  is  neceetary  it  to  multiply  each  term  In  the  residual  equation 
(Eq.  A. IS)  by  a  weighting  factor  proportional  to  the  velocity,  at  given  by  Eq.  A.3,  so  that,  In 
forming  the  coefficients  of  the  normal  equation  (Eq.  A.  14),  the  square  of  this  factor  would  enter 
each  product.  This  should  be  do  it  with  caution,  however,  since  it  will  give  very  large  weights 
to  data  near  the  beginning  of  the  range.  It  Is  also  clear  that  the  approximate  value  of  b  mutt  be 
known  in  advance  in  order  to  determine  euch  a  weighting  factor,  namely,  1  ♦  (b  Hi).1-1 

Finally,  It  wae  obeerved  that  the  value  obtained  for  a  wae  quite  sensitive  to  fluctuations  In  the 
data  near  the  end  of  the  record  of  obeervatlona.  This  Is  to  be  expected,  of  courtt,  elnce  a  rep- 
reeents  the  atymptote  of  tht  velocity.  An  unutuilly  large  deviation  near  the  end  of  the  record 
or  a  gradual  changt  In  the  epeed  of  the  camera  would  exert  an  exaggerated  influents  on  the 
value  obtained  for  a.  Whenever  a  set  of  data  encompasses  a  vary  limited  range  of  radii,  if  le 
probably  better  to  atslgn  the  velocity  of  sound  psrmansntly  to  t  and  to  determine  only  b  and  c 
from  the  data.  If  It  should  be  desired,  the  values  of  any  of  ths  parameters  can  be  fixed  In  ad¬ 
vance  and  only  ths  remaining  ones  determined  according  to  tho  leuet- square*  criterion. 
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Defense  Nuclear  Agency 

6801  Telegraph  Road 
Alexandria,  Virginia  22310-3398 


errata 


19  October  1995 


MEMORANDUM  FOR  DEFENSE  TECHNICAL  INFORMATION  CENTER 
ATTENTION:  OCD/Mr.  Bill  Bush 


SUBJECT:  Classification  Review  of  AD-363573L 


The  Defense  Nuclear  Agency  Security  Office  has  reviewed  and 
declassified  the  subject  report  (AD-363573L,  WT-613) . 

Distribution  statement  "A"  (Approved  for  Public  Release) 
applies . 

FOR  THE  DIRECTOR: 


t 


A 


JOSEPHINE  h.  WOOD 
Chief,  Technical  Support 
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tKKAIA 


Lir  4  tfc  W\> 

fee  Mi  o*\ 


Defense  Special  Weapons  Agency 
6801  Telegraph  Road 
Alexandria.  Virginia  22310-3398 


OPSSI 


JUN  I  I  1997 


MEMORANDUM  FOR  DISTRIBUTION 

SUBJECT:  Declassification  Review  of  Operation  IVY  Test 

Reports 


The  following  31  (WT)  reports  concerning  the  atmospheric 
nuclear  tests  conducted  during  Operation  IVY  in  1952  have  been 
declassified  and  cleared  for  open  publication/public  release: 

WT-602  through  WT-607,  WT-609  thru  WT-618,  WT-627  thru  WT- 
631,  WT-633 ,  WT-635 ,  WT-636,  WT-639,  WT-641  thru  WT-644,  WT-646, 
and  WT-649. 

An  additional  2  WTs  from  IVY  have  been  re-issued  with 
deletions.  They  are: 

WT-608 ,  WT-647 . 

These  reissued  documents  are  identified  with  an  "Ex"  after 
the  WT  number.  They  are  unclassified  and  approved  for  open 
publication . 

This  memorandum  supersedes  the  Defense  Nuclear  Agency,  ISTS 
memorandum  same  subject  dated  August  17,  1995  and  may  be  cited  as 
the  authority  to  declassify  copies  (of  any  of  the  reports  listed 
in  the  first  paragraph  above. 


